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Large amplitude coherent state superposition generated by a time-separated 
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Theoretical analysis is given for a two-photon subtraction from a continuous wave (cw) squeezed 
vacuum with finite time separation between two detection events. In the cw photon subtraction 
process, generated states are inevitably described by temporal multimode states. Our approach is 
based on analytical formulae that are mathematically simple and provide an intuitive understanding 
of the multimode structure of the system. We show that, in our process, the photon subtracted 
squeezed vacuum is generated in two temporal modes and one of these modes acts as an ancillary 
mode to make the other one a large amplitude coherent state superposition. 

PACS numbers: 03.67.Hk, 42.50.Dv 



I. INTRODUCTION 

Photon subtraction is a useful technique to condition- 
ally manipulate nonclassical state of light. To subtract 
photons from a traveling wave, one puts a highly trans- 
missive beam splitter (BS) and detects a reflected beam 
by photo-detectors. Selection of the event such that the 
detectors observe n photons in total approximately acts 
as annihilation operations a n to the initial state. Appli- 
cations of the photon subtraction have been proposed so 
far, e.g. generation of Schrodinger cat-like states [l|, dis- 
tillation of Gaussian entangled states [2], and loophole 
free tests of the violation of Bell's inequality [3, |4j. In 
addition, by adding displacement operations in front of 
the detectors, one would obtain higher flexibility of the 
output state synthesis [fa, [& 0] • 

Here the optical Schrodinger cat state means the su- 
perposition of distinct coherent states 

\C±) = -^(\a)±\-a)), (1) 



and an n-photon subtracted squeezed vacuum with even 
(odd) n is approximately equal to |C+) (|C_)). Recently, 
single-photon subtraction from a squeezed vacuum has 
been experimentally demonstrated with a pulsed [8( and 
cw [9j, |l0| squeezed vacuum, those correspond to the 
generation of |C_) with \a\ 2 sw 1. It should also be 
noted that an alternative way by using a photon num- 
ber state and conditional homodyne detection has been 
proposed and experimentally demonstrated with a pulsed 
two-photon state [1JJ. The conditional output corre- 
sponds to a superposition of displaced A-squeezed states 
D(±/3)S(r)\0) with \(3\ 2 w 1.2 which would become \C+) 
with \a\ 2 fa 2.6 after applying appropriate squeezing op- 
eration in P-direction. Once such a larger amplitude 
coherent state superposition becomes directly available, 
it would be an important resource for the linear optics 
quantum computation scheme [12j. 



For cw experiments, one has to take into account a 
problem of mode mismatch between a squeezed state 
and photon detection events. While coherence time of 
the cw squeezed vacuum generated from an optical para- 
metric oscillator (OPO) is given by the inverse of the 
OPO cavity bandwidth (C^ 1 ), the photon detection usu- 
ally occurs within almost instantaneous time duration 
(At <C CrT 1 )- As a consequence, the conditional output 
state appears locally around the photon detection time 
with a nontrivial mode function. General theory to sim- 
ulate such experiments has been developed in [lj, [l4| 
and the optimal mode functions are investigated in de- 
tail for the conditional generation of single- and two- 
photon states from a two-mode squeezed vacuum (non- 
degenerate opo) us m. 

In this paper, we theoretically investigate a two-photon 
subtraction from a cw squeezed vacuum generated by a 
degenerate OPO. For a pulsed scheme (or the original 
proposal in |l|), it conditionally generates an even par- 
ity superposition |C+) with \a\ 2 rs 1. Compared to the 
pulsed scheme, a distinct feature of the cw scheme is that 
the photon detection events generally occur in different 
times. We show that, surprisingly, with appropriate time 
difference A, a size of the cat state is drastically increased 
to \a\ 2 = 2.5 ~ 3 with the fidelity of F > 0.9. It is shown 
that the time-separated photon subtraction generates a 
nontrivial two-mode state distributed in two particular 
temporal modes which allows us to synthesize the out- 
put state and results a larger amplitude coherent state 
superposition in an appropriate temporal mode. Also, 
contrasted to the previous theoretical analyses of the cw 
schemes [Tj, [Ij, [lji [l6| , our approach is mostly analyt- 
ical, which is useful for intuitive understandings of the 
multimode structure in the generated state and how the 
size of the cat state is increased via quantum interfer- 
ences. 

The paper is organized as follows. In Sec. II, be- 
fore treating cw sources, we briefly review the properties 



of the photon subtracted squeezed state and coherent 
state superposition. In Sec. Ill, we discuss the cw time- 
separated photon subtraction in detail and show how we 
can extract a larger coherent state superposition from a 
temporal multimode photon subtracted state. We also 
discuss the physical insight of the same scheme from an 
alternative point of view in Sec. IV, where we show that 
the two types of quantum interferences play crucial roles 
to increase the size of the superposed coherent states. 
Section V concludes the paper. 
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II. SYNTHESIS OF APPROXIMATE 
COHERENT STATE SUPERPOSITIONS 



FIG. 1: Fidelity vs the size of the cat state a for F 
|(C"+|^2>| 2 (solid line) and F = \(C+\2PS)\ 2 (dashed line). 



Generation of a Schrodinger cat-like state by pho- 
ton subtraction from a traveling wave squeezed vacuum 
was proposed in [l|] which consists of a small reflectance 
beam splitter and, in ideal, a photon number resolv- 
ing detector at the reflected port. When the detector 
counts n photons, the transmitted state is transformed 
to \nPS) ex a n S(-r)\0), which approximates \C+) (|C_)) 



for even (odd) n, where S(r) = exp[l 



it 2)1 i s the 



squeezing operator and — r represents the squeezing in P 
direction of the phase space. Asymptotically, for large n, 
the size of a increases as \a\ 2 — n. 

Since S(rYaS(r) — acoshr — <p sinhr, n-photon sub- 
tracted state a n S(—r) |0) is rewritten as a squeezed state 
of a superposition of even or odd number states as 



iPS) = 



1 



VKi 



a n S(-r)\0) 



—=S(-r)(a cosh r + fi f sinhr) n |0) 

V-A/n 



S(-r) 



K\n) + c r n _ 2 \n-2) 

n : even, 
(c r n \n) + c r n _ 2 \n-2) 

n : odd, 



c5|o» 

clll)) 



(2) 



where Af n is a normalization factor and each of c\ is a 
function of r. For example, a single-photon subtracted 
state is exactly equivalent to a squeezed single-photon 
state. 

Meanwhile, it has recently been predicted that if one 
could arbitrarily synthesize the superposition ratio Cj in 
the last line of Eq. ([2]), approximate \C±) with larger a 
would be obtainable 0, [HI ■ To discuss it more precisely, 
let us consider the state 



the optimal coefficients Cj's maximizing F are propor- 
tional to those of the X-squeezed \C±) 

1 



S(r)\C ± 
up to n. Here 
S(r)\±a) = 



A/± 



S(r) (\a) ± | - a)) , 



(5) 




m),(6) 



where A = tanh r and H m (x) is a Hermite polynomial. 

For example, when n = 2, the two-photon subtracted 
state is given by 

\2PS) ex a 2 S(-r)\0) 

= S(-r) sinhr (V2 sinhr] 2) +coshr|0)] .(7) 

The squeezed superposition state \ip2) is on the other 
hand found from Eq. ([5]) to be 



|^ 2 > oc S(-r) 



(1-A 2 )a 2 -A 
V2 



|2> + |0) 



(8) 



The fidelities between these states and |C+) are plotted 
in Fig. ([1]) where note that, for given a, the squeezing 
parameter r is optimized to maximize the fidelities. It is 
clearly shown that one can obtain \C+) with more than 
twice larger average power (\a\ 2 = 2.5 ~ 3) if it is possible 
to synthesize the number state superposition. 

In the next section, we show that, with a cw squeezed 
vacuum source, one can synthesize the superposition by 
simply having a time separation between two photo- 
detection events. 



\ipn) = S(-r)(c n \n) +Cn_ 2 |n- 2) •• •), 



and suppose we can arbitrarily set c^'s. Since the fidelity 
between \C±) and \ip n ) is described as 



F = \(C±\S(-r)(c n \n) 



(3) III. TIME-SEPARATED PHOTON 
SUBTRACTION FROM A CW SQUEEZED 

VACUUM 

In this section, we derive an analytical expression of the 

(4) time-separated two-photon subtracted squeezed vacuum. 
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FIG. 2: (Color online) Schematic of the two-photon subtrac- 
tion from a cw squeezed vacuum. 



Schematic is shown in Fig. O A cw squeezed vacuum 
generated from an OPO with the bandwidth (o is split 
via a small reflectance BS with the reflectance R and the 
reflected fraction is guided into two photon detectors in 
path B and C through a half BS {R x = 1/2). Choosing 
the event such that two photons click each of detectors 
at time t — t\ and £2, where |^2 ~ t\\ is comparable or 
smaller than (^ > one can conditionally obtain a two- 
photon subtracted state as an output in path A, which 
is measured by a homodyne detector. 

For cw sources, it is important to take into account the 
mode mismatch between the squeezed vacuum and the 
photon detection events since the squeezed vacuum is a 
multi-mode state distributed within a time Cq~ , while 
the photon detections happen almost instantaneously 
[la , [lj, [jjl, [l6[ . The conditional output we want to ob- 
serve is therefore generated in a temporal mode localized 
around t\ and ti. To detect such state, the homodyne 
measurement might be a time integrating detection with 
an appropriate mode function filter. In other words, the 
homodyne measurement extracts one particular temporal 
mode from the cw signal. 

In the following subsections, we first give a slight mod- 
ification of the usual input-output theory of the degener- 
ate OPO, and then describe the cw photon subtraction 
process in Schrodinger picture. Since the purpose of this 
section is to describe the structure of the time-separated 
photon subtracted state, we assume that all elements in 
the scheme are lossless and there is no technical noise. 
We show that, in an appropriately filtered conditional 
state, the number state superposition is well synthesized 
via the time difference of photo-detection events. 



Squeezing operation via an ideal degenerate 
OPO 



Let us start by briefly reviewing the input-output the- 
ory of the degenerate OPO [17| . Denote the positive- 
frequency part of the field operator in a single transverse 



mode of the optical field with a continuous spectrum by 



1 f 

a(t) = — / dfla(uJo + Sl)e 
2tt ./-00 



Here a(u>o+£l) is the annihilation operator for the Fourier 
component at angular frequency luq + Q, where loq is the 
center angular frequency of the OPO. The operator obeys 
the continuum commutation relation 

[a(u> + n), a 1 (wo + n')] = 2n5 ( n - n ')- ( 10 ) 

The time dependent field operator a(t) is defined in the 
interval (—00,00), and obeys the commutation relation 

[a(t), at (0] = *(«-*'). (11) 

These operators should be moved to a rotating frame 
at the center frequency luq, 



A(t) = a(t)e iUot = -!- / </<) UO)r 



-iQt 



(12) 



where A(£l) = a(uio + ft). 

Following [17], an input-output relation of the lossless 
OPO is described by the Bogolubov transformation 

S f A A(Q,)S A = n(il)A{il) + i/(f2)it(_n), (13) 



where 



A*(«) 



Co 2 + e" 



n 2 



2Coe 



(14) 
(15) 



(Co-^) 2 -e 2 ' 

Here £0 — 7t/2 corresponds to the bandwidth of the 
OPO output, and 7t and e are the leakage rate of the 
OPO output mirror and the parameter proportional to 
the nonlinear coefficient of the OPO crystal and the 
pump amplitude, respectively. 

The above transformation is rewritten in a useful form 
by introducing the input field annihilation operator 

A e (Q) = A(ny e(n \ 



with 



tan0(O) 



2c n 



Crj 



ft 2 



(16) 



(17) 



which obeys the communication relation 

[A e (a),Al(Cl')} = 27n5(f2 - fi')- 



(18) 

The Bogolubov transformation in Eq. (|13[) is then rewrit- 
ten as 



& A A(fi)§ A = jifflMto) + p(n)ij(-n), (19) 



where 



ji{n) = n{VL)e- l6i - n) = 



v{U) = iy{Q)e~ ie ^ = 



and C± = Co ± e- 



Co 2 



2Coe 

^(a + ^xe + ff 2 ) 



,(20) 
,(21) 



B. Schrodinger picture of the two-photon 
subtraction from a cw squeezed vacuum 

Let us now turn to the photon subtraction operations. 
Denote the time dependent field operator in paths A, B, 
and C (in a rotating frame) by A(i), B(t), and C(t) , 
respectively (see Fig. [2j. A cw squeezed vacuum state 
generated from a lossless OPO is described by §a\Oa)- 
After beam splitting it via the two BSs (R and i?i), the 
quantum state distributed in three paths is described 
by VbcVabSa\Oabc) where Vab is a beam splitting 
operator transforming the field operator, for example, 
V\ B B(t)V AB = -y/RA{t) + VT^RB(t). Suppose the 
first and second photons are detected at the time t = t\ 
and t 2 in paths B and C, respectively, while the detec- 
tors project the state onto vacua in all other time. We 
assume that the detector's time resolution is instanta- 
neous i.e. enough shorter than the cavity lifetime (^ . 
After such event, the conditional output state in path A 
is projected onto 

\Pc W ) CX {0Bc\C{t2)B{ti)VBcVABS A \QABc) 

R 



where 



^ exp In y/i — R 

xi(i 2 )i(i 1 )^ A |0 A ) 



dtA\t)A{t) 



(22) 



where note that R\ = 1/2. In the limit of 
small R, the exponential term is approximated to be 
exp[ln VI - R J dtA^ (t)A(t)] ~ 1 and thus the state is 
described by a two-photon annihilated squeezed vacuum 



\p cw ) <x A(t2)A(ti)S A \0 A ) 



(23) 



By use of the property of the squeezing operation in 
Eq. (fi"9|) . we further obtain 

\ Pcw ) ex A{t2)A{ti)SA\0 A ) 

= S A \ dtv{t - t 2 )A\{t) / dtv{t - ti)A\{t) 



+ dtfL(t-t 2 )P(t-ti))\0A), 



(24) 



where 



-tn(t-to) 



(25) 



i r°° 

/i(i-io) = 77- / dfi/2(fi)£ 
27r J-oo 
1 /-oo 

*(t-to) = TT" / dm?(n)e~ in(t ~ to) . (26) 

27T ./_„ 



Defining the normalized temporal function 

1 






&i(±. 1 

2 U- Ch 
the output state is simply described as 



(27) 
(28) 



\p cw ) ex S A (M v A\A[ + Fa) \°a), (29) 



A 



F A 



t - 



dti>{t-ti)A\{t), (30) 



dtp,(t-t 2 )v(t-h), (31) 



and A = \t 2 — ti\ is the time difference between the pho- 
ton detection events. Note that, in the limit of small 
e, the temporal mode function is approximated to be 
i)j{t) rs v / Co e_l ' '*' which is a reasonable approximation 
in the realistic experiments [9|, [lCJ. Now, since xl)(t — ti)'s 
are nonorthogonal to each other, it is useful to introduce 
an orthonormal basis in these two-dimensional tempo- 
ral modes. Let us here recall that we will observe only 
a single-temporal mode at the final homodyne detection 
step, which should be a highly nontrivial state. One of 
the natural choices of the such mode is to include the first 
term in the right hand side of Eq. (j2"5|) as much as possi- 
ble. Then we arrive at the basis consisting of symmetric 
and anti-symmetric orthonormal functions defined by 



*±(*) 



j^t - t 2 ) ± jj(t - h) 
a/2(1±/a) ! 



(32) 



where 



I a = / dtip{t-t 2 )ip{t-ti) (33) 

J — OO 

and \E l +(t) is expected to be the final mode, i.e. the filter 
function of the homodyne detection. Corresponding field 
operators and state vectors are also defined as 



.4 



t - 



dW±{t)Al(t), 



(34) 



and 

Al\n ± ) = V^TT|n+l ± ), \0 A ) = |0+)|0_)|0j[). (35) 
Then \p cw ) is now expressed as 

2+,0-> 



S A (K(1 + I a ) 



Vat I V2 

AC(1-/a) 



0+,2_)+Fa|0+,0_)HOi), 



where 



V2 ' J+ 
S A \p+-)\0 A ), 

N = Nf^ + ID + Fl 



(36) 



(37) 



it it 



Note that, in Eq. (|36)l . the cross-term AI^-Al is vanished 
due to the bunching-like interference of these creation 



operators. We will discuss such quantum interference 
again in the next section. 

The conditional output state described by Eq. |36|) 

is understood as the squeezed state of |/j_| )|0^) where 

|/5-i ) is the superposition of the two-photon states occu- 
pying the temporal mode *ff+(t) or *&_(£), respectively, 
and the vacuum state. Note that when two photons are 
detected at the same time (A = 0), the second term of 
\p+~) vanishes (mode &-(t) cannot be defined) and thus 
the ideal two-photon subtracted state is localized into 
the temporal mode ^f + (i). Moreover, as will be shown, 
when A increases, one can observe a larger amplitude 
coherent state superposition in mode \£+(t). In the fol- 
lowing, therefore, we choose ^+(t) as the filter function 
for the temporal mode of the LO field in the homodyne 
detector, i.e. the detector's observable is described by 
the quadrature operator 



X 



HD 



dt^+(t)X(t), 



(38) 



(in practice, the integration is carried within a finite time 
width V (>Co _1 ))- 

Let us see the reduced quantum state in mode >F + (£). 
We first trace out mode A from \p cw ). The reduced state 
is then given by 



P+- 



= Tr A [\p cw )(p cw \] 



(39) 



where 5+ and <S_ are the non-unitary Gaussian opera- 
tions acting on modes \I/+ and , F_ separately. Each of 
S± can be regarded as a single-mode squeezing opera- 
tion although the process includes a small coupling with 
thermal environment (for derivation and characteristics 
of S± , see Appendix A) . 

Let us further trace out mode *B-(i) from /5_| We 

then obtain 



P+ 



Tr_[p+_] 
C^+\<t>){<f>\+Cj+\0){0\, 



where 



c 2 |2 H 
1 



co|0 H 



Ia), 



y/KH I V2 



-FA|0^ 



and 



k p = ^(i+i A r+Fi 



(40) 



(41) 



(42) 



As shown in these equations, the output quantum state 
observed by the homodyne detector with the filter func- 
tion ^f+(t) is the statistical mixture of the squeezed \<j>) 
and vacuum. The ratio of the statistical mixing is given 
by 



C,i, 



A^ 



±AC 2 (1 + Ia) 2 



n 



A/*(l 



-iD + K 



(43) 
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FIG. 3: The upper and lower lines show Cp and C v , respec- 
tively with |e|/£o=0.1 (solid lines), 0.3 (dashed lines), and 0.5 
(dotted lines). 
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FIG. 4: The upper and lower lines (at (oA = 0) show cq 
and c 2 , respectively with |e|/£o=0.1 (solid lines), 0.3 (dashed 
lines), and 0.5 (dotted lines). 



and C v = 1 — C,/,. 

In Figs. [3] and 2J we plot C& 



On 



and 



that 



characterize p + , as functions of A. For large delay (A 3> 
I/Co) j the output state goes to the mixture of two- and 
zero-photon state as 



P+ 



1 



S + |2)(2| + |0)(0| 



(44) 



which was already pointed out in the scheme of generat- 
ing multi-photon state from a cw non-degenerate OPO 
source Q- 

In our scheme, on the other hand, the intermediate 
delay region (A ~ (q 1 ) is particularly interesting. Since 
C v is almost negligible as shown in Fig. [31 p + can be 
approximated to be 



p+« s + \</>)(</>\. 



(45) 



On the other hand, Fig. [4] clearly shows that the tuning 
of A around (q 1 provides us a wide controllability of the 
superposition ratio of two- and zero-photon states in \(f>). 
As discussed in Sec. II, this allows us to generate larger 
coherent state superposition. 

To see the superposition property in phase space, we 
have calculated the Wigner function of p + without any 
approximation (see Appendix B). Typical Wigner func- 
tions are shown in Figs. 02a)- (c) for different A. These 
figures clearly show that the increase of A induces a 
larger size superposition. Note that the fidelity between 
the state in Fig. E^b) and \C+) with \a\ 2 = 2.6 is 0.946. 
Fidelities between p + and |C+)'s with different ampli- 
tudes are shown in Fig. [5] as a function of A. We see 
that our time-separated two-photon subtraction tech- 
nique allows us to generate a larger superposition such as 
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FIG. 5: (Color online) Wigner functions of p + . |e|/Co = 0.27 
and (a) C A = 0, (b) CoA = 1.4, (c) CoA = 2.4. The fidelity 
between (b) and the ideal cat state with \a\ 2 = 2.6 is 0.946. 




FIG. 6: Fidelity between p+ and the ideal cat state with 
H 2 =2, 2.5, 3, and 3.5 (from left to right). |e|/Co = 0.27. 



\a\ 2 > 2.5 with high fidelity (0.9 >), that is not possible 
with a normal single-mode (pulsed) two-photon subtrac- 
tion. 



IV. QUANTUM INTERFERENCE OF THE 
SUBTRACTED PHOTONS 



In the previous section, we have observed a bosonic 
bunching-like effect between A' + and A_ that originates 
from the two-photon subtraction (see Eq. (|36|) ). In this 
section, we again discuss the cw two-photon subtraction 
from the viewpoint of quantum interferences. It is natu- 
ral to expect that the state synthesis is achieved by some 
quantum interference. Generally, quantum interference 
occurs when the quantum process or states are going to 
be intrinsically indistinguishable. We concretely show 
that, in our case, two types of such interferences con- 



tribute to the size increase of the cat state. 

Let us reexamine the cw state \p C w) = 
-T=A(t2)A{ti)S A \0) and consider from which mode 

the subtracted photons {Afo) and A(t±)) originates, 
&+(t) or *£_ (t). This is clarified by defining the following 
annihilation operators 



A t± = 



7f( i(i2) 



±A(h) 



(46) 



where A t+ (A t _) is the annihilation operator which sub- 
tracts one photon from the squeezed vacuum in mode 

This is confirmed by directly substituting them as 






UN- Al]SA]0) 

a + S A h + )\0-)\0 A )-cr-S A \0 + )h-)\Q A ), 



(47) 



where we clearly see that A t+ and A t _ change the states 
only in \&+(£) and ^-(t), respectively. Moreover, from 
the first to second line, we observe a bunching-like quan- 
tum interference, i.e. the possibility such that each pho- 
ton is subtracted from each of i &±(t) is vanished. Note 
that 



|7±) 



1_( r fl±e~(- A 1±e-<+ A 



"-. I V C- c 



and 



a 2 ± = 2 



C- 



l±e^- A 



l±e^- A 

C- 



<4 



(48) 



l±e" 


C+A\ 


l±e 


/ 
-C+A 



<4 



(49) 



The final state to be measured is obtained by tracing 
modes A and *_(*) out. If |0_) and |7_) in Eq. (g7|) are 
orthogonal (distinguishable) , the final state is reduced to 
be a statistical mixture of |0+) and |7+). These terms 
are, however, nonorthogonal, i.e. partially indistinguish- 
able. One therefore may expect that some of the quan- 
tum coherence between |0+) and I7+) remains. This is 
the second type of quantum interference. To see this 
clearer, let us first simplify the equation by taking the 
approximation of the small squeezing limit. Neglecting 
the higher order of e/£o, we have 



\ Pcw ) « a + S A ^A{ 2 + l)\i) ..;,j0 :-\i)y- 

- <T _S A (^-At*+\ j l);.|l). 0,, ( r,i)) 



where 



1 ± e-« oA 



S± 



= 1± 



CoAe~ c ° A 



(51) 
(52) 



Note that S + is around 1 ~ 1.25 and in practice, e/Co 
is often regarded as an effective pumping parameter nor- 
malized by the OPO threshold. We should argue that 
the state where the two photons are subtracted at the 
same time t\ (A = 0) is also given by 



1 



A 2 (h)s\o} «s 



Co' 



-AV 



i 



10) 



(53) 



which corresponds to the original photon subtraction 
scheme proposed in |l[ and here we call it a time- 
degenerate two-photon subtracted state. 

Let us trace modes A and $f-(t) out and see the re- 
duced state with finite A. Approximating the squeezing 
operation acting on modes \l/ + (£) and \£_(t) to be a sepa- 
rable unitary squeezing operator S+ <X> S— (see Appendix 
A), we obtain the reduced state as 



time of the squeezed vacuum, the conditional output is 
appeared within two temporal modes, as a squeezed su- 
perposition of 1 20), 1 02), and |00). Such superposition 
and a careful choice of the single mode function 1 &+(t) 
allow us to synthesize the even photon number superposi- 
tion of the conditional output state and with appropriate 
parameters it results a generation of cat-like states which 
have more than 90% fidelity with the coherent state su- 
perposition of \a\ 2 > 2.5. We have also discussed the 
same issue from the viewpoint of quantum interference, 
which reveals how the conditional output state of the 
time-separated photon subtraction is deviated from that 
of the time-degenerate one, due to the quantum inter- 
ferences. Our theoretical approach provides analytical 
expressions of the states, which would be further useful 
to investigate an intuitive physical picture of more com- 
plicated multi-mode cw quantum states. 



APPENDIX A: SQUEEZING OPERATION ON 

TEMPORAL MODES *±(t) VIA AN OPTICAL 

PARAMETRIC OSCILLATOR 



p+n(l-C v )\$)($\ + C v S + \0)(0\S. 



t 

+ ■ 



(54) 



where C v = ^(1 - I A ) 2 / {g(l + I 2 A ) + e~ 2 ^ A } and 



|<I>> x,S'_ ( ^±i f + 2 + l j |0) - l.ml) 



CoA 
2 



S+|0>,(55) 



these correspond to Eqs. (|40p and (|4"Tj) . Therefore, the 
component of the cat-like sate |$) produced by the time- 
separated two-photon subtraction stems from the quan- 
tum interference (superposition) of the time-degenerate 
two-photon subtracted state and a squeezed vacuum. 
That is, the engineering of the state indistinguishability 
in mode \P_(£) allows us to synthesize the superposition 
ratio in mode \l/+(i), in other words, provides the size 
controllability of |$) in mode *f>+(t). Note that, on the 
other hand, |0_) and J7_) are only partially indistinguish- 
able. Their distinguishable part creates the statistical 
mixture term S+|0)(0|S+ with the factor of C v . 



In this Appendix, we discuss the input-output relation 
of the OPO from the input state prepared in modes \&± (t) 
to the output state in modes \P± (i) . Namely, we look at 
the completely positive trace preserving map 



S p+- = Tr A S A (A 



l0iX0.il) S A 



(Al) 



In the following, we assume that the OPO is lossless 
and we often use the mode functions ^-t^) that are 
the Fourier transformed expressions of \l/± (t) . 

Let us define the complete orthogonal set 
{y±(Q),y> { j v \Q)} i in frequency domain where Vj"\t) 
corresponds to the mode function for the vacuum input. 
Applying it into the annihilation operator A(Q), we have 

A(Q) = *+(fi)i+ + *_ (0)i_ + J2 *i V) ( fi )4 (t,) , (A2) 



CONCLUSION 



where 



In conclusion, we have theoretically investigated the 
time-separated two-photon subtraction from a continu- 
ous wave squeezed vacuum. In a single-mode theory, a 
two-photon subtracted squeezed vacuum is regarded as 
a squeezed state of d^\2) + Cq\0) in which, however, the 
superposition coefficients c!^ and Cq are not optimal to 
maximize its fidelity to a coherent state superposition. 

In case of a cw scheme, one needs a multi-mode theory 
and we showed that when the time difference of the two 
photo-detection events is finite but within the coherence 



A^ = - j ,/<ni/! ( (j).i ( ()). 



(A3) 



and again the superscript (v) means that its initial state 
is a vacuum. 



From the Bogolubov transformation of the OPO 
A out (fl) = S A A in {n)S A , which is defined in Eq. 1[I9"]). 
we can describe the input-output relation of the OPO 
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FIG. 7: Squeezing parameters |r+| (solid lines) and |r_| 
(dashed lines). From the lower to upper curves, |e|/Co = 0.1, 
0.3, 0.5, respectively. 
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FIG. 8: Thermal photons n+ (solid lines) and n_ (dashed 
lines). From the lower to upper curves, |e|/Co = 0.1, 0.3, 0.5, 
respectively. 



with respect to A± as 



1 



2 out 

1 

2: 



— OO 



(A4) 



drw* ± (n)A out {n) 

i f°° / 

— / dn\^ ± {n)\ 2 f/2(fi)i± + p(fi)i± f 

Note that we have used the relation 

dQj2(Cl)^X(fl)^ _(0) 

= — / dOP(fi)#+(O)#_(Q) = 0, (A5) 

which implies that the OPO does not couple the modes 
V^+(f2) and ^_(J7). It therefore means that the map S 
can be decomposed as 



1 
2^ 



5 = «S 4 



(A6) 



The concrete expressions of the maps S± are given as 
follows^ Since the OPO includes only Gaussian opera- 
tions, S± can be fully characterized by the real matrices 
S± , Y± , which describe the input-output relation of the 
covariance matrix as 

r5 rt = J s^r ln s , ± + y±, (A7) 

where Tin and r± ut are the covariance matrices of the 
input and output states, respectively, consisting of the 
variances of the quadratures, e.g. 



r, 



{xl 



(-Xin-Pin) 
(PI) 



From Eq. (|A4|) . we find 



(xr 2 ) = Gi 2 (x? a )+t( F ±- G f^ 



T x 



(R 



out 2\ 



r .±2/p2\ 
U P Vin/ 



li^-cr) 



( ^outpou t) = G ± G ±(X ln P in ), 

where 

1 r°° 
Gi = — J dn 9 {n)\^ ± {n)\ 2 

4C 2 6 2 ] 



(A8) 

(A9) 

(A10) 
(All) 



M/(1±/a)2tt 



dn 



1 ± cos Aft 



1 /■ 

g^ = — y dfiflf-^nji^cn)! 

4C 2 T 2 ] 



(C+ + n 2 ) 1/2 (CH. + ^ 2 ) 3/2 ' 

(A12) 

2 



AC(1±/a)2tt 



dn 



1 ± cos An 



(Cl + n 2 ) 3 / 2 (c 2 _ + o 2 )V2' 



F 



A' 



— / d0 5 2 (0)|* ± (0)p 



4«e 



2,2 



dn- 



1 ± cos AO 



C 2 e 2 {i±(i + C-A) e -c- A } 



Fi 



1 

2^ 



C1AC(1±/a) 

drj ff - 2 (^)|* ± (fi)| 5 



(A13) 



2/2 



4Co^ 



1 



rffi 



1 ± cos AO 



AC(1±/a)27t7_ 
C 2 e 2 {l±(l + C + A)e-^ A } 



and 

Then we obtain 



(C^ + tt 2 ) 2 ^ 2 +n 2 )- 



i a + » 2 

C 2 +« 2 ' 



5± = 

*± = 



G x ° 
G± 

p± r<± 2 






(A14) 

(A15) 

(A16) 
(A17) 



The above OPO process includes squeezing operation and 
the coupling with a squeezed thermal environment. The 
effective squeezing parameters r± can be defined as 



where 



r± 



1, G x 
2 Gt 



(A18) 



and how the process is deviated from a unitary squeezing 
is roughly estimated by looking at the numbers of the 
thermal photons n± 



n± = G X G% - 1. 



(A19) 



Figs. [7] and [8] plot r± and n±, respectively, which clearly 
show that the OPO process with respect to each of ^f± 
can almost be regarded as a single-mode unitary squeez- 
ing. 



APPENDIX B: WIGNER FUNCTION OF p+ 



The Wigner function of p+ in Eq. (|40j) . where S+ is 
applied on the non-Gaussian state C ( f ) \(j>)((j)\ + G„|0)(0|, is 
derived from its characteristic function. The character- 
istic function X+(w, v) is calculable with the help of the 
Bogolubov transformation in Eq. (|A4[) as 



X<t,{u,v) = \l- —7x{u 



Xo(u,v) 



x exp 
exp 



v i ) + i( U 2 + vr 



i(u 2 + v 2 ) 



;(u 2 + v") 



(B2) 
(B3) 



X+(u,v) = Tr 


P e i{uX + +vP+) 






= Tr 


'\p cw }(Pc W \e i(uJC++VP+) ~ 




= ( P+ -\(0 A \ (& A e««*++" p +)§ A ) \P+-)\0 A ) 


= G X0 (G x u, G p v ) 
xXo(V^-G+ 2 




u,^F+-G+ 2 v\ 


+CvXo (Jf+u, ^F+v) , (Bl) 

I 



Then its Fourier transformation gives the Wigner func- 
tion 



W+(x,p) 



n^jF x F+ 

Gp ly/2. 
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C2 C 



- V2( 
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2 I U X x 4 
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(2p 2 -F+) 



exp 
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Fi 



OCoCz i 



V2c 2 c + V2c 2 )- ^2L_ (2a * F+) 



F 



x 



„2,^+2 



\j v Up 



-2^X 



+ 2n+2 



F + F 



(2x 2 -F+)(2p 2 -F+; 



(B4) 
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